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Abstract: We study the thickness of the confining flux tube generated by a pair of sources 
in higher representations of the gauge group. Using a simple geometric picture we argue 
that the area of the cross-section of the flux tube, as measured by a Wilson loop probe, 
grows logarithmically with source separation, as a consequence of the quantum fluctuations 
of the underlying k-string. The slope of the logarithm turns out to be universal, i.e. it is the 
same for all the representations and all the gauge theories. We check these predictions in 
a 3D Z4 lattice gauge model by comparing the broadening of the 1-string and the 2-string. 
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1. Introduction 

In the physics of quark confinement there are many indications that the gauge field responds 
to a static q source separated from a conjugate q source by a large distance L by forming 
a colour fiux tube which behaves as a string with energy V{L) ~ L, where the string 
tension depends on the representation ^ of the quantum numbers carried by the source. 

If the gauge group is SU(iV) there are infinitely many irreducible representations at our 
disposal to cast the sources. However, for large separations, no matter what representation 
is chosen, depends only on the A''— ality k of i.e. on the number (modulo A^) of 
copies of the fundamental representation needed to build ^ by tensor product, the reason 
being that all representations with the same k can be transformed into each other by the 
emission of a proper number of soft gluons. As a consequence the heavier strings decay into 
the string of smallest string tension. The corresponding string is referred to as a k-string. 

The spectrum of k-string tensions has been extensively studied in recent years, in the 
continuum as well as on the lattice [^-[0]. In this paper we want to explore another 

facet of k-string physics, related to the quantum fiuctuations of these objects. 

It is widespread belief that the fiux tube generated by a pair of sources q, q in the 
fundamental representation is in the rough phase. This means, as explained long ago by 



Liischer, Miinster and Weisz |17], that the colour fiux tube broadens as the sources are 
separated. More precisely, the area of a cross-section of the fiux tube should increase 
logarithmically with separation. This phenomenon can be seen as a consequence of the 
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quantum vibrations of the underlying 1-string describing of the infrared properties of this 
flux tube. 

Since the k-strings can be seen as bound states of k 1-strings, it is interesting to ask 
if there is any fundamental obstruction for a logarithmic broadening of the k-strings. At 
first sight it is not clear whether in the IR limit the only relevant degrees of freedom are 
the transverse displacements of a single string or whether we have to take into account 
new degrees of freedom describing the possible splitting of the fc-string into its constituent 
strings. In other terms one would like to know whether the string binding energy damps 
significantly the quantum fluctuations of the free strings. 

The answer we find to the above questions is surprisingly simple: not only the logarith- 
mic broadening occurs in the k-strings for any k, but it turns also out that the coefficient 
of the logarithm does not depend on the specific representation of the source nor on its 
N-ality and is universal, i.e. it is the same in all gauge theories. We check these predictions 
in a 3D Z4 gauge model, which is the simplest environment where a non-trivial 2-string 
can live. A brief report of our work has been presented in [18|. 

Measuring the thickness of the flux tube and in particular its dependence on the 
separation of the sources is very challenging from a computational point of view. In SU(A^) 
gauge theories the error bars are too large to draw deflnite conclusions [19|. As a matter of 
fact, an uncontroversial observation of logarithmic broadening has been only made in 3D 
Z2 gauge model, thanks to the efficiency of the Monte Carlo algorithms for its dual, the 
Ising model |20|. New results on the thickness of the SU(2) flux tube near the deconflning 
point exploiting the integrability of the underlying 2D Ising model appeared recently pl| . 
The numerical data on the logarithmic growth of the mean squared width of the flux tube 
associated to the Z4 strings presented here constitute a new important support on the 
expected quantum behaviour of the flux tube. 



2. The rough phase 

In the strong coupling phase of whatever gauge theory in three or four space-time dimen- 
sions the flux tube joining a quark pair has a constant width for large inter-quark distances. 
As the coupling constant decreases, the flux tube can undergo a roughening transition. The 
rough phase is characterised by strong fluctuations of the collective coordinates describing 
the position of the underlying string and the mean squared width of the flux tube diverges 
logarithmically when the inter-quark distance goes to infinity [^]. 

The square width of the flux tube generated by a planar Wilson loop Wf{C) in the 
fundamental representation is defined as the sum of the mean square deviations of the 
transverse coordinates hj{S,i,(,2) (i = I7 • • • — 2) of the underlying string, i.e. 

= T E / - h^f) , (2.1) 

j=i ■'^ 

where is the planar domain bounded hy dO. = C , A = d^^ its area and h°- are the 
transverse coordinates of the equilibrium position. The vacuum expectation value is taken 
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with respect to the two-dimensional field theory describing the dynamics of the flux tube. 
It is widely believed that the roughening transition belongs to the Kosterlitz-Thouless 



universality class |22|. As a consequence, it is expected that this fleld theory flows, in the 



infrared (IR) limit, toward the massless free- field theory described by the Gaussian action 

D-2 „ 

^^« = f E / E d,h,d^h,. (2.2) 

j=l -'^ Ai=1.2 

In such a limit the mean square width can be easily evaluated in terms of the free Green 
functions Gn{C,0 = iKOHC)) as 

= ^ ^ d^C 1^ {GniC, a - Gni^, C + e)) , (2.3) 

where e is a UV cut-off. The action ( |2.2D is conformally invariant, hence the integration of 
the finite part G{^, ^') cannot depend on the size of the domain Q but only on its shape. 
The logarithmic growth of w'^ comes from the UV divergent part. This can be simply 



understood as follows [20|. The conformal invariance of the theory implies the scaling 
property 

Gn{^,0 = GnAs^,s^'), (2-4) 

where s > is an arbitrary real number and 0,s is the scaled domain. On the other hand 
in the UV limit — > the Green function diverges logarithmically 

Gn{tC + e) = --^\oge + ... . (2.5) 



|2.4| ) and (2^) agree only if the cut-off appears exclusively in the ratio e/R, where i? is a 



typical linear size of the domain. Adding this piece of information to (|2.3| ) yields 

w^ = ^^logiR/Rfm, (2.6) 

where the UV cut-off has been absorbed in the scale Rf{i^), thus the absolute value of 
this physical scale cannot be determined by this conformal approach. On the contrary the 
ratios of these scales for different domains are calculable functions of the shapes pO| , p3| . 

It is clear from this approach that the generalisation to the k-string crucially depends 
on its IR limit. If, for instance, it fiuctuated as a single string it would suffice to put 
(T (Tk in ( |2.2| ) and modify the other formulae consequently. There are other possibilities, 
however. We shall see in the next section that in Z) = 3 dimensions there is a geometric 



approach, first advocated in |17], which can be unambiguously extended to the k-strings. 



3. K-strings and minimal surfaces 

To inspect the width of the flux tube generated, in a 3D gauge system, by a static source in 
the fundamental representation one can start, following jl^], from the connected correlator 

p.,^ {Wf{G)Wf{c))-{Wf{G)){Wj{c)) 
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where Wf{C) and Wf{c) are Wilson operators for the loops C and c, which are concentric 
circles of radii R and r < R lying on parallel planes at a distance h, as shown in Fig. |l|. 
The above quantity can be seen as a measure of the density of the colour flux generated 
by the source C and felt by the probe c. 

To study the flux generated by sources in a generic 
representation we will take W,^{C) instead of Wf{C), 
but keep the probe always in the fundamental. 

The mean square width of the flux tube is defined 

by 

2 _ fh'^Po^{h)dh 



(3.2) 






Figure 1: Parallel circular loops 



The insight of ||l^ was to observe that the quantity 
Pf{h) can be described in the effective string picture in 
terms of the world-sheet of a Nambu-Goto string con- 
necting the two circles. This defines a typical Plateau 
problem of minimal surfaces, which can be evaluated in 
the infrared limit by a saddle-point approximation. More specifically, one can write 

( Wf{C) Wf{c) ) - { Wf{C)){Wf{c) ) oc eM-^MR, r, h)) , (3.3) 

where A{R, r, h) denotes the area of the connected minimal surface having C and c as 
boundaries. 

The choice of the Nambu-Goto (NG) ac- 
tion for the effective string is the simplest one, 
however there are well known problems for its 
quantisation procedure. This theory is fully 
consistent only at critical space-time dimen- 
sions (D=26)^, but there is strong evidence 
that the first few terms in the \/{aA) ex- 
pansion of the NG action are universal up to 
the order O (^r^p- (see |2^, and references 
quoted therein, for a detailed discussion on 
this point). In particular, the 0(1) term co- 
incides with ( p.2[ ), thus it is reasonable to ex- 
pect that in the infrared limit R ^ oo with 
fixed r and h the result should not depend on the specific choice of NG action. Note how- 
ever that the size r of the probe cannot be too small, otherwise the string picture would 
not be valid. 

The minimal surface connecting the two circles is a surface of revolution about the 
symmetry axis (here chosen to coincide with the x axis). If we denote by y{x) the y > 
section of the surface with the {x,y) plane, the area is given by 

A 




Figure 2: Catenoid surface 



2-K 



yVy 



r,2 



1 dx . 



(3.4) 



^Polchinski and Strominger developed an effective string tfieory wfiicfi avoids quantisation problems 
in non-critical physical dimensions. 
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The variational condition = yields the equation 1 + y'^ = yij. The general 
solution is ^ 

y(x) = — cosh.uj(x — xq) , (3-5) 

that is, the surface of revolution is a catenoid (see Fig. ^). The integration constants must 
obey 

ii =—coshu;xo , r = — cosh. uj{h — xq) . (3-6) 
UJ u> 

This results in the following expression for the minimal area: 



moreover Eq. (^^) allows to express /i as a decreasing function of u> 

, , rarccosh(i?a;) — arccoshfro;)] 
h(uj) = ^ ^ ^ — ^ , 3.8 

a; 

thus we can regard u; as an integration variable for actually computing nP' . Using the 
trivial inequality coshx > 1 (Vx) we get a minimal value for a;, that is a maximal allowed 
value for /i: ^ 

^min — ~ ^ hqjiax — h{\/r) . (3-9) 

r 

Now we can write explicitly, for the mean square width, 

2 iT K^? exp[-aA[R,r,h{uj))] \h'{uj)\duj 

J'l'exp[—aA{R^r,h{oj)y\\h'{uj)\doj ^ ^ 

r 

This quantity approaches a logarithmic curve for large R. Indeed, this can be seen by 
using the asymptotic expansion 

arccosh(cjr) ^ 21og(u;r) {uor oo) , (3-11) 



and inserting it into ( |3.8D , which yields the approximate solution found in [17|: 

a; ~ i log(i?/r) . (3.12) 

The condition ur ^ 1 becomes 

log{R/r) :^ h/r , (3.13) 

always fulfilled in the large R limit. Note that r cannot be too small at fixed R and h. In 
this limit, a Gaussian distribution is found for the transversal density, whose width grows 
logarithmically with R: 



p(«)(/i) 



oc exp 



\og{R/r) 



aw^ = ^log{R/r), (3.14) 



which is almost identical to (|2.q ). Here the UV cut-off is replaced by the size r of the 
probe. 
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The result found here can be generahsed to Wilson loop W^{C) in any representation 
(where we are no more guaranteed that the infrared Gaussian free string limit is valid). 
The world-sheet of the k-string can be seen as some bound state of k 1-string world-sheets. 
In this case a local minimum solution exists for sure, in which one of these fundamental 
sheets gives rise to the catenoid with the probe Wf{c), while the other k — 1 just lie flat 
on the loop surface. We then have immediately 



(fc) 



>(K) 



(1) 



(3.15) 



where the stability of the k-string implies that the exponent is always negative. The 
resulting width, then, appears to be exactly the same as for the fundamental representation. 



3.1 Special configurations 

In some special cases, for a very narrow range of h, there is, beside the above general 
solution, another minimal surface made by a first catenoid composed by the k-string world- 
sheet. 

At a suitable distance d it splits into a disk orthogonal 
to the symmetry axis made by the world-sheet of the {k — 1)- 
string and a second catenoid in the fundamental representa- 
tion which reaches the probe c. The position of the intermedi- 
ate disk is not arbitrary, but it is dynamically determined by 
the balance of the tensions at the string junction. In the most 
general configuration depicted in Fig. ^ a A;2-string propagat- 
ing along a catenoid decays into a /ca-string forming a disk and 
a /ci-string forming another catenoid. The angles of the string 
junction are given by 




Figure 3: The balance of 
the string tensions 



cos 6*1 



a 



k2 



(3.16) 



and cyclic permutations of the indices for the other angles. 

Implementing these kinematic constraints on the parameters a;o,xo and uji,xi of the 
two catenoids yields 



ujQ R = cosh ujoXq , ujq p = cosh ujq{xq — d) = 1/ sin 6i , 
Wi r = cosh (xi + d — h) , uji p = cosh loiXi = 1/ sin 02 . 



(3.17) 
(3.18) 



where d is the distance of the disk from the source and p its radius. We can invert these 
relations and express d and h as functions of p 



d{p) = 
h{p) = p sin 



p sin^i arccosh(ii//9 sin^i) — /9 sin^i arccosh(l/sin0i) , (3.19) 
62 arccosh(l/ sin02) — P sin02 arccosh(r//9 sin 6*2) -|- d{p) . (3.20) 



Both d{p) and h{p) turn out to be increasing functions of p. Since h > d we have p > r 
which implies 62 > -f, as Fig. ^ shows, then ( 3.16| ) tells us that this special surface is 
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permitted only if the kinematic constraint 



4, > <yl, + < (3-21) 

holds, i.e. only when the binding energy of the /c-string is sufficiently small. On the other 
hand, the fact that the argument of arccosh must be larger or equal to 1 fixes the maximum 
Pmax = T I sin Q2 ■ Thus the minimal surface exists only in the range 

h{r) < h{p) < h{r/ sin 62) , (3.22) 

which turns out to be in general very narrow. These special configurations should produce 



a spike in the distribution of the flux density at the distance fixed by (3.22). 

The hight of such a peak is determined by the total area of the surface. Unfortunately 
this area depends significantly on r, hence it is rather unclear how to extract physical 
information on the flux tube, which should not depend on the size of the probe. 

In the N = 4:, k = 2 case, which corresponds to our simulations, this special solution 



should obey the constraint ( 3.21 ) i.e. a2 > V^cr. In our numerical data taken in that 



region we did not find signs of such a singular behaviour. 
4. Z4 gauge theory and its dual 

The laboratory where we study the physical properties of the 2-string is a Z4 gauge model 
in three dimensions, defined by the standard plaquette action on a cubic lattice. With Z4 
as gauge group, there exist two k-strings in the system: the fundamental string and the 
2-string (related to the double- fundamental representation / ® /). Moreover, since there 
are no more representations than one in the same N-ality class, the system does not exhibit 
any meta-stable string that could spoil the results at finite R. 

It is well known that this model, as any three-dimensional abelian gauge model, admits 
a spin model as its dual (see for instance ) and any physical property of the gauge system 
can be translated into a corresponding property of its spin dual. From a computational 
point of view it is of course much more convenient to work directly on the spin model where 
powerful non-local cluster algorithms can be applied. 

In our case the dual is a spin model with global Z4 symmetry which can be written as 
a symmetric Ashkin- Teller (AT) model [p7[], i. e. two coupled, ferromagnetic, Ising models 
defined by the two-parameter action 

Sat = - ^ [/?(o-a'0-j/ + r^Ty) + a{ar,ayTr,Ty) , (4.1) 

where ax and Tx are the Ising variables (cr, r = ±1) associated to the site x and the sum is 
over all the links (xy) of the dual cubic lattice. The phase diagram of this model has been 
studied long ago |2^, p9| . 

The global Z4 symmetry of the action is generated by the transformation 

a ^ -t; t ^ a . (4.2) 
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The model has also an independent Z2 symmetry generated by 



(7 ^ r , (4.3) 

related to the charge conjugation of the corresponding dual model. 

It is customary to represent the Z4 symmetry through the multiplication table of the 
fourth roots of the identity (^j = ±i,±l). It would be a simple exercise to rewrite the AT 
model with this kind of variables introducing the complex field 

^^ = e-'^{a^ + iT^), (4.4) 

where the phase factor is chosen such that the Z2 symmetry defined above becomes the 
complex conjugation * -^-^^ ^f*. The Z4 symmetry, instead, corresponds to the multiplication 
by i. 

The standard application of the duality transformation would lead to a gauge model 
on the dual lattice with a {.^j}-valued field on the links and with an action containing 
plaquette operators and their squares: 

S = -J2 h^^Up) + C2iUpf] . (4.5) 
p 

We think it is interesting (and more purposeful for our work) to show that the AT model 

is dually equivalent to two coupled Z2 gauge models. 

The starting point is to rewrite the Boltzmann factors associated to the links in the 
known form, namely the Z2 character expansion 



.pT.cTy ^ ^^^^^^ ePf^M^a^ayYi-y) (4.6) 



g/3r.r, ^ /Sillh2^ ^ ^'i3u^xy)^^^^^y(xy) (4.7) 



V^xy)=±l 



^aaxayTxTy ^ ^ ^ p^^y) ^^^^^r^^yfi^y) (4.3) 

P{xy)=±^ 

with 7= — log v'tanh7 and ^ = 5 — 2 • 

Now the sum over the site variables a and r can be explicitly performed, yielding for 
each node x the two conservation laws 

Yl ii^i^y) + Pi^y)) = 2 , (z>(^y) + p{^y)) = mod 2 (4.9) 
y y 

or, equivalently, in the multiplicative form 

Y[P{yx)P{yx) = ^Yi^{yx)P{yx) = '^- (4-10) 
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The canonical partition function becomes 

Zat oc ^ expy^ 



/3 {l^{xy) + T^{xy))+ " P{xy) 



(4.11) 



The apex in the sum over configurations indicates that they must obey the constraints 



( 4.10 ). A way to solve them in an infinite lattice is suggested by the usual duality transfor- 
mation of the three-dimensional Ising model: it is sufficient to consider the two composed 
link variables lJ'[yx)P{xy) ^-^id V{xy)P{xy) the plaquette variables of the dual lattice. More 
precisely we solve the above constraints with the Ansatz 

P(xy) = ^(xy)Up , V{xy) = ^{xy)yP ; P{xy) = ^(xy)UpVp , (4-12) 

here e is an arbitrary sign variable (e = ±1), P is the plaquette dual to the link {xy). The 
plaquette variables are defined through the product of their boundary links, namely, 

Up= Y,Ue,Vp= Y,Ve, (4.13) 

e&dP eedP 

where i are links of the dual lattice, of course. It is a straightforward exercise to verify 
that the above Ansatz solves identically Eq.s ( 4.10| ). The sum over the e variables is 



unconstrained and can be performed at once, leading to 

Zat oc Yl n [PiUp + Vp)+ a UpVp] . (4.14) 

{Ue=±l,Ve=±l} P 

As the last step, we can rewrite this partition function in the usual Boltzmann form by 
defining 

^^b{Up+Vp)+aUpVp ^ ^Qgj^ []3{Up + Vp)+ a UpVp] . (4.15) 
where a, b, c are suitable coefficients. Solving for a, b, c we get 

= cosh(2 P + a) cosh(2 P - a) cosh^ a (4.16) 
e46 ^ cosh(2 j + a) 



cosh(/3 — a) 
4a _ cosh(2 f3 + a) cosh(2 f3 — a) 



cosh'^ a 



(4.18) 



therefore we can recast Z^t as the partition function of two coupled Z2 gauge systems 

Zat(/9,«)oc Y1 J]e''(^^+^^)+"^^^^, (4.19) 

{Ui=±iye=±l} P 

where the duality transformation & : (a, /?) (a, b) from the AT couplings to the gauge 
couplings can be explicitly written, using ( 4.171 ) and ( [4.18| ), as 



Ij / (coth /3 -1- tanh/3 tanh a)(coth /? -|- tanh/?coth a) \ f4 20) 

4 V 2 + tanh a + coth a ) 

1, / 1 -|- tanh^ /5 tanh a \ 

--In ^ . (4.21) 

4 V tanh^ /? tanh a / 
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As a non-trivial check of these formulae one can verify that the ^ transformation is an 
involutory automorphism, i. e. = 1, as required for any duality transformation. 



4.1 Dual of a Wilson loop 

The duality transformation maps any physical observable of the gauge theory into a cor- 
responding observable of the spin model. In particular it is well known that the Wilson 
loops are related to suitable twists of the couplings of the spin model. More specifically, 
let us consider a Z^r spin model in 3D. Let ^ = e^'^'^l^ be the generator of this group. A 
/c-twist of the link (xy) in the spin action is defined by the substitution /? ^(i only in 
the selected link. Denoting with Zi^r^y^u the spin partition function modified in this way, 
one can easily prove the identity 

{Up^^ ) gauge = Z(^r^y-^^k/Z , (4.22) 

where P is the plaquette dual to {xy) and C/jf'' is the plaquette variable in the irreducible 
representation of Z^r characterised by the integer /c = l,2,...,A^. A simple check is the 
following. The twist of the spin variable associated to a single node is equivalent to twisting 
all the links incident to this node; on the other hand the twist of a single spin can be re- 
absorbed in the invariance of the measure of the partition function. On the gauge side, 
this corresponds to the fact that the product of the plaquette variables lying on the six 
faces of the cube dual to the selected node is identically equal to 1. 

Repeating the above construction for a suitable set of plaquettes we can construct in 
this way any Wilson loop or Polyakov-Polyakov correlator in any representation and its 
map into the spin model. 

We want to fit this procedure to the AT model. In this case twisting a link corresponds 
to associating to it an anti-ferromagnetic coupling, i. e. (3 ^ —(3. The character expansion 
(|4.6| ) for an anti-ferromagnetic Boltzmann factor becomes simply 

,-P^.^y = ^^l ^ /x^.,)e^^(-) (a.a,)'^<-) . (4.23) 

Notice that in the AT model there are three different couplings associated to a single link. 
Which of them do we have to twist in order to build on the gauge side the plaquette in 
the fundamental (i. e. A; = 1) representation? To answer this question it suffices to twist, 
for instance, the variable ax — > —^x- From the point of view of the symmetries of the AT 
model, such a twist corresponds to a Z4 generator (^]^) followed by a charge conjugation 
( |4.3| ); this corresponds to ^x ™- terms of the complex variable defined in ([4.4|), 

hence it is associated to the fundamental representation. On the other hand this change of 
sign yields the twist of two couplings for each link incident to x: the quadratic a coupling 
and the quartic coupling. Therefore the plaquette in the fundamental representation is 
simply obtained by multiplying the Boltzmann factor by Ai(xy> P{xy) = where we used 
(4.12). In this way, the identity outlined above can be written explicitly as 



( Vp ) gauge = { e-2(/3+--^-.)---. (4.24) 



-10- 



generalising the known dual identity of the Ising model. Similarly, flipping the signs of 
both spins Ux and Tx we get the plaquette variable in the k = 2 representation as ( UpVp ). 
Combining together a suitable set of plaquettes we may build up any Wilson loop or 
Polyakov-Polyakov correlator with /c = 1 or A; = 2. 

The reason why we insist in writing this model in terms of Ising variables is that 
here one can easily apply a very efficient non-local cluster method |31] which generalises 



in a straightforward way the one commonly used in Ising systems, based on the Fortuin- 
Kasteleyn (FK) cluster representation. Moreover, it has been built a very powerful method 
to estimate Wilson loops ( ) based on the linking properties of the FK clusters for 
each FK configuration generated by the above-mentioned algorithm one looks for paths in 
the cluster linked with the loop 7. If there is no path of this kind we put = 1, otherwise 
we set = 0. This method leads to an estimate of ( ) with reduced variance with 
respect to the conventional numerical estimates. 

5. Monte Carlo simulations 

We performed a Monte Carlo analysis on the AT model, for both the fundamental and 
the double-fundamental strings, at the (confining) coupling (a,/3) = (0.0070,0.1975) (for 
which we have measured the string tensions a? a = 0.01560(1) and a? 02 = 0.0210(6) pO[|). 

On the lattice, the quantity P^{h) is measured by placing a square Rx R loop W{R) 
on a plane in the desired representation and taking the probe as a plaquette operator, 
parallel to the loop and lying on its axis at a distance h. The actual single measurement 
took into average also the four planar neighbours of the plaquette in the central position, 
in order to enhance the signal; this operation does not spoil the results since we dealt with 
large values of R. 

According to the recipe for embedding the presence of a Wilson loop (in the represen- 
tation ^) directly into the action as a series of frustrated links, the algorithm has only to 
measure the expectation value ( Up )^-'^ of the probe plaquette in the fundamental rep- 
resentation, where the superscript outside the average symbol denotes the fact that the 



action is modified to include the loop W{C) according to (3.1) 



P^{h) oc ( C/^ constant . (5.1) 

In terms of AT variables, measuring a plaquette operator translates to measuring the 
coupling energy on its dual link (xy) (the sign variable e is equal to -|-1 for every link except 
those dual to the loop, where its value reflects the applied frustration) 

( )Yauge 0^ ( ecFxCTy ) AT + Constant . (5.2) 

The algorithm used for the analysis uses a cluster update method basically similar 
to the standard Fortuin-Kasteleyn cluster technique: each update step is composed by an 
update of the a variables using the current values of the r as a background (thus locally 
changing the coupling from /3 to /3 it a according to the value of r^Ty on the link {xy)), 
followed by an update of the r's using the a values as background. 
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Note that since the large loop is automatically handled by the update procedure and 
the probe is a loop of side 1 there is no need to implement expensive topological analysis on 
the configuration to get the measures, allowing us to reach a performance, with the lattices 
we used, of about 0.35 seconds per update/measurement step (as on a single Intel ® Xeon 
3.2 GHz 64-bit processor). We used a cubic lattice with side L = 80 (we found there are 
no finite size deviations up to i? ~ L/2) and measured times the plaquette operator as 
discussed. Configuration results have been then packed in groups of 25 in a binning fashion 
to estimate variances. We performed the measurements on loop sides R = 11, 13, ... ,41 
with a statistics of 220875 measures (with independent configurations for each R) for the 
fundamental representation and 470275 for the k = 2 case. 



6. Results 

The argument based on the properties of the minimal surfaces describing the world-sheet 
of the underlying confining string suggests a logarithmic growth of the mean square width 
of the flux tube for both the fundamental and the double-fundamental {k = 2) string. 
Also the numerical values of the two functions w^{R) at fixed R should coincide, but this 
prediction is easily spoiled by artifacts of the short-range difference between a bound state 
of two strings and a fundamental string. 
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Figure 4: Typical {k — 2, R — 23) transversal flux density distribution, as measured with the 



quantity in Eq. (5.2). Note the mismatch with the Gaussian form and the bad background estimate 
that would follow from the Gaussian assumption for P{h) 
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Figure 5: Flux tube width, k ~ 1 versus k — 2 



We found that the measured values of P<%{h) do not fit too well to a normal distribution 
as approximately expected and lead to bad width estimates (see Fig. so we used the 
numerically integrated quantities instead, as in eq. (3.2). We performed the integration 
by carefully choosing a cutoff value hmax ■ since the background value must be subtracted 
from the transverse flux density functions, the results are very sensitive to the choice of a 
cutoff. We took hmax = 12. To estimate this background value, we looked for a plateau in 
the region \h\ > h'^j^^^ separately for each R, and took as final value the broadest averaged 
result in the range fi^-^ = 18, . . . , 24. 

By fitting the functions 2(^)5 R — Rmin with Rmin an appropriate distance 
cutting off non-IR contributions, to the functional form 



aw^{R) 



27r 



log R + c 



(6.1) 



(see Fig.^) we found that the fundamental string width, as well as the 2-string, show the 
expected logarithmic growth with the appropriate universal multiplicative factor (reduced 
were, respectively, 1.22 and 2.68), but the value of c differs measurably in the two cases: 
ci = 0.4002(20), C2 = 0.4512(11); this discrepancy is probably due to some interaction 
between the fundamental world-sheets in the k = 2 case. 

Our work reinforces the numerical evidences of the predicted logarithmic broadening 
of the flux tube width [20|, extending them with high precision to the case of Z4, the 



simplest gauge group with more than one kind of k-string. In particular, the main result is 
that, when dealing with a non-fundamental string, its effective width still grows with the 
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logarithmic law. This suggests that here the assumption of a free massless string behaviour 
for large R holds from the k = 1 case. 



7. Conclusions 



In the context of SU(A^) gauge theories we have studied the thickness of the confining 
flux tube generated by a pair of static sources in higher representations as probed by a 
Wilson loop in the fundamental representation, whose extension is small compared with 
the source separation. Generalising a simple confining string picture proposed long time 
ago by Liischer Miinster and Weisz ||l^ we argued that mean square width w"^, when 
measured in terms of fundamental string tension units a, grows logarithmically with the 
source separation R in a manner which is universal, i.e. 



D 



a w 



2tt 



log R/Rc . 



(7.1) 



The reference scale Rc cannot be directly determined by the underlying string model and 
it is the only place where one can envisage a dependence on the source representation. 

We performed a careful verification of these predictions in the case of a 3D Z4 gauge 
theory, which is the simplest gauge system where a 2-string forms. 

To reach the required high precision of the numerical data, the simulations where 
actually performed in the dual version of the system, which turns out to be a symmetric 



Ashkin- Teller model |27]. This choice allowed us to use efficient non-local cluster algorithms 
1 31]. The results of this analysis compare very favourably with the above predictions. 
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